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Abstract:- In this paper, the DE optimization technique is effectively used to solve the optimal power 
flow problem by incorporating Facts device i.e. SSSC to enhance the performance of the power system. 
The standard IEEE 30-bus test system is considered to examine proposed approach without and with 
SSSC FACTS device. Results show that proposed DE algorithm gives better solution than other 
algorithms to enhance the system performance with SSSC device. 
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I. INTRODUCTION 

In today's highly complex and interconnected power systems, there is a great need to improve electric 
power utilization while still maintaining reliability and security. While power flows in some of the transmission 
lines are well below their normal limits, other lines are overloaded, which has an overall effect on deteriorating 
voltage profiles and decreasing system stability and security. Because of all that, it becomes more important to 
control the power flow along the transmission lines to meet the needs of power transfer. On the other hand, the 
fast development of solid-state technology has introduced a series of power electronic devices that made FACTS 
a promising pattern of future power systems. Power flow is a function of transmission line impedance, the 
magnitude of the sending end and receiving end voltages and the phase angle between voltages. By controlling 
one or a combination of the power flow arrangements, it is possible to control the active as well as the reactive 
power flow in the transmission line [1]. With FACTS technology [2], such as Static Var Compensators (SVCs), 
Static Synchronous Compensators (STATCOMs), Static Synchronous Series Compensators (SSSCs) and 
Unified Power Flow Controller (UPFC) etc., bus voltages, line impedances and phase angles in the power 
system can be regulated rapidly and flexibly. Thus, FACTS can facilitate the power flow control, enhance the 
power transfer capability, decrease the generation cost, and improve the security and stability of the power 
system. 

In this paper, SSSC FACTS controller are incorporated to solve an optimization problem with different 
objectives such as minimization of cost of generation, real power loss, voltage profile enhancement and 
improvement of voltage stability L-index as these are the basis for improved system performance. The 
Differential Evolution (DE) algorithm is used effectively to solve the optimal power flow problem, it results 
great characteristics and capability of determining global optima, by incorporating a set of constraints including 
voltage stability and FACTS device. In order to calculate the power losses and check the system operating 
constraints such as voltage profile, a load flow model is used. An existing Newton-Raphson load flow algorithm 
is introduced [2]. This model is further modified to incorporate SSSC FACTS device into the network and DE 
technique is applied to the modified model to enhance the performance of the power system. Thus, effectiveness 
of the proposed method was tested on standard IEEE 30-bus test system and comparison was made on the 
performance of system with other OPF methods. 

The organization of this paper is as follows. Section 2 addresses the Computation of Voltage Stability 
Index (L-index). FACTS controller is explained in Section 3. Mathematical formulation of optimal power flow 
problem is given in section 4. Differential Evolution Algorithm Optimization Process is represented in section 5. 
The overall computational procedure is given in the section 6. The simulation results on test system are 
illustrated in section 7. Finally, the conclusion is given in section8. 

This document is a template. An electronic copy can be downloaded from the conference website. For 
questions on paper guidelines, please contact the publications committee as indicated on the website. 
Information about final paper submission is available from the website. 
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II. Voltage Stability Index (L-index) Computation 

The voltage stability L-index is a good voltage stability indicator with its value change between zero 
(no load) and one (voltage collapse) [3]. Moreover, it can be used as a quantitative measure to estimate the 
voltage stability margin against the operating point. For a given system operating condition, using the load flow 

(state estimation) results, the voltage stability L-index is computed as [4]: 



L r 



8 V 

U ■' Vj 



j = g + l,...,n (1) 



All the terms within the sigma on the RHS of equation (1) are complex quantities. The values of F /7 are 
btained from the network F-bus matrix. 

For stability, the index L must not be violated (maximum limit=l) for any of the nodes j. Hence, the 
global indicator L describing the stability of the complete subsystem is given by maximum of Lj for all j 
(load buses). An L -index value away from 1 and close to indicates an improved system security. The 
advantage of this L -index lies in the simplicity of the numerical calculation and expressiveness of the results. 

III. FACTS CONTROLLERS 

FACTS controllers are able to change in a fast and effective way, the network parameters in order to 
achieve better system performance. FACTS controllers [5,6] such as phase shifter, shunt, or series compensation 
and the most recent developed converter-based power electronic controllers, make it possible to control circuit 
impedance, voltage angle and power flow for optimal operation performance of power systems, facilitate the 
development of competitive electric energy markets, stimulate the unbundling the power generation from 
transmission and mandate open access to transmission services, etc. The benefit brought about by FACTS 
includes improvement of system behavior and enhancement of system reliability. However, their main function 
is to control power flows. 
3.1. Static Synchronous Series Compensator (SSSC): 

A SSSC [7] usually consists of a coupling transformer, an inverter and a capacitor. The SSSC is series 
connected with a transmission line through the coupling transformer. 

It is assumed here that the transmission line is series connected via the SSSC bus j. The active and 
reactive power flows of the SSSC branch i-j entering the bus j are equal to the sending end active and reactive 
power flows of the transmission line, respectively. In principle, the SSSC can generate and insert a series 
voltage, which can be regulated to change the impedance (more precisely reactance) of the transmission line. In 
this way, the power flow of the transmission line or the voltage of the bus, which the SSSC is connected with, 
can be controlled. 
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Fig. 1: Equivalent Circuit of SSSC 

The equivalent circuit of SSSC is as shown in the Fig.l. From the equivalent circuit the power flow 
constraints of the SSSC can be given as: 

p.. =vf £„. -Vyj(g, cos6l +b l] m$,)-vy„(£ t cos(6j. -0J+b t sin(6>. -0J) ( 2 ) 

Q„ =-v, 2 b„ -Ws* sjn 3, -h C0f 4)- WS,y *HO t -0 K )-b t cos0 ( -6J) ( 3 ) 
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where g v + fr v =UZ„,g u =g il ,b u =b ll ,gj l =g v ,bj l =b t (6) 

The active and reactive power flow constraints is: 
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where r-- and {J t are specified active and reactive power flows. 



specified 
i 

specified -. s~\ specified 

and Q j. 
The equivalent voltage injection V se Z.0se bound constraints are as : 

x / min < y <ymax 

se se se x ' 

mm <0 <6> max (10) 

se se se \ ±yj s 

IV. MATHEMATICAL FORMULATION OF OPF PROBLEM 

Mathematically, the OPF problem with FACTS is solved to minimize fuel cost of generation 
maintaining thermal and voltage constraints can be formulated as follows [14]-[22]: 

NG 

Minimize F = (£(a i P* i +b i P Gi + C i ) (11) 

1=1 
The minimization problem is subjected to following equality and inequality constraints 

4.1 Equality Constraints: These are the sets of nonlinear power flow equations that govern the power system: 

^-^-ZW|^K|cos((9 s -^+^) = (12) 

7=1 

Q Gl -Q Dl + t J \V,\^ j \\Y lj \sm(e ij -S,+S J ) = (13) 

7=1 

where P Gi and Q Gi are the real and reactive power outputs injected at bus I , the load demand at the same bus 
is represented by P Dj and Q Di , and elements of the bus admittance matrix are represented by \Y r and ■■ . 

4.2 Inequality Constraints: These are the set of constraints that represent the system operational and security 
limits like the bounds on the following: 

1) generators real and reactive power outputs: 

P Gi ^ P ci^ P ci ,i=l..., ng (14) 

at ^Qa^Qci ,i=l-..,ng (15) 

2) voltage magnitudes at each bus in the network: 

vr^v^vrj^i-.^g (i6) 

3) transformer tap settings: 

T i mm <T i <Tr\i = l,...,NT (17) 

4) reactive power injections due to capacitor banks: 

Q™r<Q a <Q™r,i = l...,CS (18) 

5) transmission lines loading: 

S.<S. max ,*' = l,...,n/ (19) 

6) voltage stability index: 

Lj i <Lj™i = l,...,NL (20) 

7) SSSC device constraints: 

SSSC Series voltage source magnitude 

V mm <V <V max (21) 

se se se v y 

Series voltage source angle 

er<e se <er (22) 

The equality constraints are satisfied by running the power flow program. The generator bus real power 
generations ( P ' • ), generator terminal voltages ( V ■ ), transformer tap settings ( T t ), the reactive power 

generation of capacitor bank ( Q ci ), P.. and Q-. of SSSC are control variables and they are self-restricted by 

the representation itself. The active power generation at the slack bus ( P ), load bus voltages ( V u ) and 
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reactive power generation ( Q t ), line flows ( S t ), and voltage stability (L ■ )-index are state variables which are 
restricted through penalty function approach. 

V. DIFFERENTIAL EVOLUTION ALGORITHM OPTIMIZATION PROCESS 

A differential evolution algorithm (DEA) is an evolutionary computation method that was originally 
introduced by Storn and Price in 1995. DEA uses rather greedy selection and less stochastic approach to solve 
optimisation problems than other classical EAs. There are also a number of significant advantages when using 
DEA. 

A Initialization: 

In the first step of the DEA optimization process, the population of candidate solutions must be 

initialized. Typically, each decision parameter in every vector of the initial population is assigned a randomly 

chosen value from within its corresponding feasible bounds: 

„ . f 7nax _ mitiL^ 

X j j., nTmL i^ /iJ .\X j Xj ) (23) 

where i — l,...,NP and j = \,...,D. X;} is the initial value (G=0) of the jth parameter of the i' 
individual vector, xf"" 1 and xf am are the lower and upper bounds of the f decision parameter, respectively. Once 
every vector of the population has been initialized, its corresponding fitness value is calculated and stored for 
future reference. 

B. Mutation: 

The DEA optimisation process is carried out by applying the following three basic genetic operations; 
mutation, recombination (also known as crossover) and selection. After the population is initialised, the 
operators of mutation, crossover and selection create the population of the next generation P (G+1) by using the 
current population P {G) . At every generation G, each vector in the population has to serve once as a target vector 
Xf G \ the parameter vector has index i, and is compared with a mutant vector. The mutation operator generates 
mutant vectors (V/ G) ) by perturbing a randomly selected vector (X rl ) with the difference of two other randomly 
selected vectors (X r2 and X r3 ). 
V&* = x?? +- F X^ - X%* , i = L„, N p (24) 

Vector indices rl, rl and r3 are randomly chosen, which rl, rl and rl {l,...,NP} and rl ^ r2^ rl ^ i. Xrl, Xrl 
and Xr3 are selected anew for each parent vector. F is a user-defined constant known as the "scaling mutation 
factor", which is typically chosen from within the range [0, 1+] . 

C. Crossover: 

In this step, crossover operation is applied in DEA because it helps to increase the diversity among the 
mutant parameter vectors. At the generation G, the crossover operation creates trial vectors ({/,•) by mixing the 
parameters of the mutant vectors ( V,) with the target vectors (Xj) according to a selected probability distribution: 
•?P if rand j (0,1) < CR or j = s 

X;/ otherwise / 05 n 

The crossover constant CR is a user-defined value (known as the "crossover probability"), which is 
usually selected from within the range [0, 1]. The crossover constant controls the diversity of the population and 
aids the algorithm to escape from local optima, rand, is a uniformly distributed random number within the range 
(0,1) generated anew for each value of j. s is the trial parameter with randomly chosen index {1,...,D}, which 
ensures that the trial vector gets at least one parameter from the mutant vector. 

5.4 Selection: Finally, the selection operator is applied in the last stage of the DEA procedure. The selection 
operator chooses the vectors that are going to compose the population in the next generation. This operator 
compares the fitness of the trial vector and the corresponding target vector and selects the one that provides the 
best solution. The fitter of the two vectors is then allowed to advance into the next generation according to 
equation (26): 



fl? C 



1 X^ otherwise 



(26) 

The DEA optimization process is repeated across generations to improve the fitness of individuals. The 
overall optimization process is stopped whenever maximum number of generations is reached or other 
predetermined convergence criterion is satisfied. 
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VI. OVERALL COMPUTATIONAL PROCEDURE FOR SOLVING THE PROBLEM 

The implementation steps of the proposed DE based algorithm can be written as follows; 
Step 1: Input the system data for load flow analysis 
Step 2: Select FACTS device and its location in the system 
Step 3: At the generation Gen =0; set the simulation parameters of DE and randomly initialize k individuals 

within respective limits and save them in the archive. 
Step 4: For each individual in the archive, run power flow under the selected network contingency to determine 

load bus voltages, angles, load bus voltage stability indices, generator reactive power outputs and 

calculate line power flows. 
Step 5: Evaluate the penalty functions 

Step 6: Evaluate the objective function values and the corresponding fitness values for each individual. 
Step 7: Find the new generation individuals and store them. 
Step 8: Increase the generation counter Gen = Gen+1. 
Step 9: Apply the DE operators to generate new k individuals 
Step 10: For each new individual in the archive, run power flow to determine load bus voltages, angles, load bus 

voltage stability indices, generator reactive power outputs and calculate line power flows. 
Step 11: Evaluate the penalty functions 

Step 12: Evaluate the objective function values and the corresponding fitness values for each new individual. 
Step 13: Apply the selection operator of DE and update the individuals. 
Step 14: Update the new generation and store them. 

Step 15: If one of stopping criterion have not been met, repeat steps 4-14. Else go to stop 16 
Step 16: Print the results 

VII. SIMULATION RESULTS 

The proposed DE algorithm is employed to solve optimal power flow problem by incorporating SSSC 
FACTS device for enhancement of system performance on standard IEEE 30-bus test system. The DE 
parameters used for the simulation are summarized in Table I. 



Table I: Optimal Parameter Settings for DE 



S.No. 


Parameters of Differential evolution 


Parameter 


values 


1. 

2. 
3. 
4. 


Population size 
Number of iterations 
Scaling mutation factor, F 
Crossover Factor, CR 


50 
250 
0.5 
0.9 



The network and load data for this system is taken from [22]. To test the ability of the proposed DE 
algorithm one objective function is considered that is minimization of cost of generation. In order to show the 
affect of power flow control capability of the FACTS device in proposed DE OPF algorithm, two sub case 
studies are carried out on the standard IEEE 30-bus system. 
Case (a): power system normal operation (without FACTS devices installation), 
Case (b): one SSSC device is installed in line connected between buses 9 and 10 with real and reactive power 

flows ( Pfi , and Q ■■ ) as + 1.25 times of base case values. The ratings of SSSC are: V se is in the range [0.001, 

0.2], 9 se is in the range [0,2 n ]. 

The first case is the normal operation of network without using any FACTS device, in second case 
optimal location of device has been considered. 

From the Table II, it can be seen that details of the control variables and the installation of SSSC in the 
network gives the best performance of the system in the network in terms of reduction in cost of generation, 
power loss reduction, maximum of voltage stability indices. It also gives that DE algorithm is able to enhance 
the system performance while maintaining all control variables and reactive power outputs within their limits. 



Table II: Optimal settings of control variables for IEEE 30-bus system 



Control 
Variables 


Limits(p.u) 


DE 
Without 
FACTS 


DE With 
FACTS 
device 

SSSC 


Min 


Max 
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Poi 


0.50 


2.000 


1.7714 


1.7615 


Pq2 


0.20 


0.800 


0.4869 


0.4911 


P G3 


0.10 


0.350 


0.2104 


0.1431 


Pg4 


0.10 


0.300 


0.1183 


0.1362 


P G5 


0.15 


0.500 


0.2127 


0.2523 


Pg6 


0.12 


0.400 


0.1200 


0.1200 


V G i 


0.9 


1.10 


1.083 


1.0736 


v G2 


0.9 


1.10 


1.0641 


1.0601 


V G3 


0.9 


1.10 


1.0365 


1.0355 


V G4 


0.9 


1.10 


1.0124 


0.9987 


V G5 


0.9 


1.10 


1.0339 


1.0319 


V G6 


0.9 


1.10 


1.0444 


1.0344 


Tap-1 


0.9 


1.1 


1.0471 


1.0317 


Tap - 2 


0.9 


1.1 


0.9167 


0.9009 


Tap - 3 


0.9 


1.1 


0.9529 


0.9666 


Tap -4 


0.9 


1.1 


0.9504 


0.9444 


Qcio 


0.0 


0.10 


0.0467 


0.0000 


Qci2 


0.0 


0.10 


0.1000 


0.0373 


QC15 


0.0 


0.10 


0.0786 


0.0394 


Qci7 


0.0 


0.10 


0.0795 


0.0549 


Qc20 


0.0 


0.10 


0.1000 


0.0076 


Qc21 


0.0 


0.10 


0.0546 


0.0988 


Qc23 


0.0 


0.10 


0.0160 


0.0381 


Qc24 


0.0 


0.10 


0.0270 


0.0525 


Qc29 


0.0 


0.10 


0.0233 


0.0463 


Cost ($/h) 




798.86 


795.7573 


Ploss (p.u.) 




0.0857 


0.0702 


Ljmax 




0.1281 


0.1271 



The convergence characteristic of the cost of generation DE without and with SSSC at optimal location 
is shown in Fig. 2. 
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Fig. 2: Convergence of cost of generation without and with sssc using de for ieee 30-bus system 




Fig. 3: L-index without and with SSSC device using DE for IEEE 30-bus system 

The Figures 4-6 show the percentage MVA loading of the lines, voltage profiles and voltage angles 
indices of buses without and with SSSC at optimal location. 
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Fig.4:Percentage MVA line loadings of IEEE30-bus system after optimization without and with SSSC using DE 
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Fig. 5: Voltage profiles of IEEE 30-bus system after optimization without and with SSSC using DE 
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Fig. 6: Voltage angles of IEEE 30-bus system after optimization without and with SSSC using DE 

A Comparison of fuel cost of generation without FACTS devices: 

The comparison of fuel cost of the proposed method with those of the methods reported in the literature 
is given in Table III. It can be seen that DE algorithm gives less cost of generation compared with the cost of 
generation obtained with other OPF methods. 

Table III: Comparison of fuel costs for IEEE 30-bus system 



Method 


Fuel Cost ($/hr) 


EP [16] 


802.907 


TS [16] 


802.502 


TS/SA [16] 


802.788 


ITS [16] 


804.556 


IEP [16] 


802.465 


SADE_ALM [17] 


802.404 


OPFPSO [18] 


802.410 


MDE-OPF [19] 


802.376 


Genetic Algorithm ($/hr) [20] 


803.050 


Gradient method [2 1 ] 


802.430 


PSO (proposed) without FACTS 


800.867 


PSO (proposed) with SSSC 


797.187 


DE(proposed) without FACTS 


798.860 


DE(proposed) with SSSC 


795.757 



VIII. CONCLUSIONS 

This paper has presented an OPF model incorporating FACTS controller SSSC using DE algorithm for 
enhancement of system performance. This model is able to solve power networks of any size and converges 
with any number of iterations and independent of initial conditions. The standard IEEE 30-bus system has been 
used to demonstrate the proposed method over a wide range of power flow variations in the transmission system. 
The results shows that proposed OPF with Static Series Synchronous Compensator (SSSC) scheme using DE is 
very effective compared to other methods in improving the security of the power system. 
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